HOMOGENIZATION OF FIRST ORDER EQUATIONS 
WITH u/e-PERIODIC HAMILTONIAN: 
RATE OF CONVERGENCE AS e -> 
AND NUMERICAL APPROXIMATION OF THE EFFECTIVE 

HAMILTONIAN. 
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, Abstract. We consider homogenization problems for first order Hamilton- Jacobi equations 

^ with u e /e periodic dependence, recently introduced by C. Imbert and R. Monneau, and also 

studied by G. Barles: this unusual dependence leads to a nonstandard cell problems. We study 
the rate of convergence of the solution to the solution of the homogenized problem when the 
parameter e tends to 0. We obtain the same rates as those obtained by I. Capuzzo Dolcetta 
and H. Ishii for the more usual homogenization problems without the dependence in u £ /e. In a 
second part, we study Eulerian schemes for the approximation of the cell problems. We prove 
that when the grid steps tend to zero, the approximation of the effective Hamiltonian converges 
to the effective Hamiltonian. 
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1. Introduction 



We consider homogenization problems for first order Hamilton- Jacobi equations with u e /e 
■ periodic dependence, namely 

^ nn / ul + H(±,Z,^,Du*) = 0, (t,x)e(0,+oc)xR N , 

K ' ' \ u e (0,x) =u (x), x£R N 

with the following assumptions on the Hamiltonian H: 
(HI) Periodicity: for any (t,x,u,p) GlxR^xRxl 1 

' H(t + l,x + k, u + l,p) = H(t, x, u,p) for any k £ Z N ; 

(H2) Regularity: H : H x R N x R x R N -> R is Lipschitz continuous and there exists a constant 
Ci > such that, for almost every (t, x, u,p) £ R x R N x R x R N 

\D {ttX) H(t,x,u,p)\ <Ci(l + H), \D u H(t,x,u,p)\ <Ci, \D p H(t,x,u,p)\ <C X ; 

(H3) H(t,x,u,p) —> +oo as \p\ —> +oo uniformly for (t,x,u) G R x R N x R; 
(H4) There exists a constant C such that for almost every (t, x, u, p) G R x R N x R x R N 

\D p H(t, x,u,p) • p — H(t, x, u,p)\ < C. 

Problem (jl.ip with H independent of t was introduced by Imbert and Monneau [11] as a 
simplified model for dislocation dynamics in material science. The complete model is introduced 
in [12] and leads to nonlocal first order equations of the type 

<+ (c(^)+M e (^)^ \Du e \ + H (^-,Du* \ =0 

where M e is a nonlocal jump operator and c is a periodic velocity. In the latter model, the level 
sets of the solution u e describe dislocations. 

Going back to (jl.ip . it was proved in [11] that, with H independent of t, 

• under assumptions (HI) and (H2), there exists a unique bounded continuous viscosity 
solution of (jl.ip : 
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under assumptions (H1)-(H3), the limit u° of u € as e — > exists and it is the unique 
bounded continuous solution of the homogenized problem 

u t ° + H(Du°) = 0, (t,x) G (0, +oo) x R N , 
u°(0,x) =u (x), xeR N , 

where the effective Hamiltonian H is uniquely defined by the long time behavior of the 
solution of 

„x / A = u t + F(x, -At + + + £>?;), (t, x) G (0, +oo) x R-^, 

1 j \ u(0,x) =0, xeR N . 

More precisely, we have the following theorem 

Theorem 1.1 (Imbert-Monneau, [H]). Let H be independent oft. Assume (Hl)-(HS) and 
uq G W l,oc (R N ). Then, as e — > 0, the sequence u € converges locally uniformly in (0, +oo) x R N 
to the solution u° of (jl.2p . where, for any p G R N H(p) is defined as the unique number A for 
which there exists a bounded continuous viscosity solution of (|1.3p . Moreover H : R N — > R is 
continuous and satisfies the coercivity property 

H(p) — * +oo as \p\ — > +oo. 

The proof in [11] is rather involved: it uses a twisted perturbed test function for a higher 
dimensional problem posed in R x R N x R. 

Under the additional assumption (H4), an easier proof of Theorem 11.11 was given by Barles, [3], 
as a byproduct of a general result on the homogenization of Hamilton-Jacobi equations with 
non-coercive Hamiltonians. 

Remark 1.2. The hypothesis (H4) which was not used in [llj guarantees the existence of a 
function such that 

HooOt, x,u,p) = lim sH(t,x,u,s~ 1 p). 
Moreover satisfies (H1)-(H3). 

In [3], thanks to assumption (H4), the equation for u e is interpreted as an equation for the 
motion of a graph: indeed, following [3j, for t G R, (x,y) G R^" 1 " 1 , (p x ,Py) G R^" 1 " 1 , let us 
introduce the non-coercive Hamiltonian F defined by 

1 ' n '' V ' Px ' Pv) ~\H 0O (t,x,y,p x ), otherwise. 

The function U e (t, x, y) := u e (t, x) — y satisfies 

' U[ + F^,^^,D x U",D y U^ =0, (t,x,y) G (0, +oo) xR N +\ 
U*{0,x,y)=u {x)-y, (x,y) eR N+1 . 

In [3] Barles proves that the sequence U e converges to the solution U° of the following problem 

U? + F{D X U°, D y U°) = 0, (t, x, y) G (0, +oo) x R^ 1 , 
U°(0,x,y)=u (x)-y, (x,y)eR N+1 , 
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where for (p x ,Py) G R Ar+1 , F(p x ,p y ) is the unique number A for which the cell problem 



(1.7) V t + F(t,x,y,p x + D x V,p y + D y V) = A in R x 



>A r +l 



admits bounded sub and supersolutions. This result makes it possible to solve the homogeniza- 
tion problem for (jl.ip : 
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Theorem 1.3 (Barles, [3]). Assume (H1)-(H4). Then the sequence u e converges locally uni- 
formly in (0, +oo) x to the solution u° of (|1.2[) . The function H(p) in (jl,2p can be cha- 
racterized as follows: H{p) = F(p,—1), where, for any (p x ,Py) € F(p x ,p y ) is the unique 
number A for which the equation (jl.7p admits bounded sub and supers olutions in M. : 



An important step in the proof of Theorem 11.31 consists of homogenizing the non-coercive 
level-set equation satisfied by l{[/ e > }. 
In this paper, we tackle two questions: 

• Is it possible to estimate the rate of convergence of u e to vP when e — ► 0? 

• Is is possible to approximate numerically the effective Hamiltonian? 

The first question was answered by Capuzzo Dolcetta and Ishii, [4] for a more classical homoge- 
nization problem: the estimate \\u e — < Ce~s was obtained for Hamilton- Jacobi equations 

of the type 

u e + H (x, -,u e ) = 0, 



where (x,y,p) — ► H(x,y,p) is a coercive Hamiltonian, uniformly Lipschitz continuous for \p\ 
bounded and periodic with respect to y; moreover, if H(x,y,p) does not depend on x, then the 
convergence is linear in e. We will show that in the present case, it is possible to obtain the 
same rates of convergence as e — > by adapting the proof in [lj using the arguments contained 
in [3] . Our main result on this topic is Theorem 12.11 in § [2j The main idea is to approximate U e 
(with an error smaller than e) by a discontinuous function U e which takes integer values where 
U e has noninteger values and which is a discontinuous viscosity solution of 

UI + F^-^DJJ^DyU^j =0, (t,x,y) e (0, +oo) x M 7V+1 . 

The latter equation has to be compared with (jl.5p . This approximation U e is obtained as the 
limit as 5 — > of 4>$(U e ) where (<j>s)s is a sequence of increasing functions. The method of 
Capuzzo Dolcetta and Ishii [1] can then be applied to U € . 
The second question was studied in [1] for equation 

u e + H (~,« e ) = 0, 

where (y,p) — > H(y,p) is a coercive Hamiltonian, uniformly Lipschitz continuous for \p\ bounded 
and periodic with respect to y; in this article, a complete numerical method for solving the 
homogenized problem was studied, including as a main step the approximation of the effective 
Hamiltonian by solving discrete cell problems. Error estimates were proved. Here, we will study 
the approximation of the cell problem (jl.7p by Eulerian schemes in the discrete torus. We have 
prefered to study the approximation of the noncoercive N + 2 dimensional problem (jl.7p rather 
than that of the coercive N + 1 dimensional problem (jl.3p because the solution of (jl.3p may not 
be periodic. In § [3l we prove Theorem 13. 1\ the discrete analogue of the ergodicity Theorems 
in P|, i.e. that there exists a unique real number A^* such that the discrete analogue of (jl.7p 
has a solution. The arguments in the proof are the discrete counterparts of those in [3]. Then, 
we prove Proposition 13.31 which states that the discrete effective Hamiltonian converges to the 
effective Hamiltonian when the grid step of the discrete cell problem tends to zero. 
To summarize, the paper is organized as follows: Section [2] is devoted to finding estimates on 
the rate of convergence as e — ► 0. Section [3] is devoted to the numerical approximation of the 
effective Hamiltonian by Eulerian schemes. Finally, we present some numerical tests in Section 

m 
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2. An estimate on the rate of convergence when e — >■ 

This section is devoted to the estimate of the rate of the uniform convergence of the solutions 
of (jl.ip to the solution of the equation (jl.2p in term of e. 

2.1. The main result. 

Theorem 2.1. Assume (H1)-(H4) and u G W 1 ' oc (R N ). Let u e and u° be respectively the 
viscosity solutions of (jl.ip and (| 1 . 2 j) . Then there exists a constant C, independent of e G (0, 1), 
such that for any T > 

(2.1) sup \u e (t,x) -u°(t,x)\ < Ce T ek. 

[0,T]xR N 

If uo is affine then 

(2.2) sup \u e (t,x) -u°(t,x)\ < Ce. 

R+xR^ 

2.2. Preliminary results. In this section we recall some results that will be used later to 
obtain error estimates. 

The assumptions (H1)-(H4) on H guarantee that F satisfies 
(Fl) Periodicity: for any (t, x,y,p x ,p y ) eRx R N+1 x R N+1 

F(t + l,x + k,y + l,p x ,Py) = F(t,x,y,p x ,p y ) for any k G Z N ; 

(F2) Regularity: F : Rx R^"^ 1 x M iV+1 — > R is Lipschitz continuous and there exists a constant 
C\ > such that, for almost every (t, x,y,p x ,p y ) G R x R N+1 x R N+1 

\D(t,x) F (.t, x ,y,Px,P y )\ < Ci{\p x \ + \p y \), \DyF(t,x,y,p x ,p y )\ < Cx\p y \, 

(F3) Coercivity: F(t,x,y,p x ,p y ) — > +oo as — > +oo uniformly for (t,x,y) G E x M^ 1 , 
< for any > 0; 

Remark that F(£, x, y, 0, 0) = 0. This and (F2) imply that for every (t, x, y,p x ,p y ) G R x R N+1 x 

R 7V+1 

(2-3) \F{t,x,y,p x ,p y )\ < d(\ Px \ + \ Py \). 

Moreover, by construction, F satisfies the "geometrical" assumption 
(F4) For any (t,x,y,p x ,p y ) G R x R N+1 x R N+1 and any A > 0, 

F(t,x,y,Xp x ,Xp y ) = XF(t,x,y,p x ,p y ). 

Assumption (F4) guarantees that (|1.5p is invariant by any nondecreasing change U — > <f(U), see 
[5] and [ID], i.e., any function V = ip(U e ), with nondecreasing is solution of 

V t + F^,^,^,D x V,D y v) =0, (t,x,y) G (0, +oo) x R N +\ 
V(0,x,y)=<p(u (x)-y), ' (x, y) G R^ 1 . 

Finally, note that (F3) and (F4) imply the existence of a positive constant C2 such that 

(2.4) F(t,x,y,p x ,0) > C 2 \p x \ for all (t,x,y,p x ) eRx M JV+1 x R^. 

In [3], in order to construct sub and supersolutions of (|1.7j> . Barles introduces for a > the 
auxiliary equation 

(2.5) W t a + F(t,x,y, Px + D x W a , Py + D y W a ) + aW a = 0, (t, x, y) G R x R^ 1 , 

with F defined by (jl.4p . and shows that if (H1)-(H4) hold true, then (|2.5p admits a unique 
continuous periodic viscosity solution. Moreover the limit of aW a (t, x, y) as a — > + does not 
depend on (t,x,y) and the half-relaxed limits of VF" — minW a provide a bounded subsolution 
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and a bounded supersolution of (jl.7p . with A = — lim Q ,_ +0 + aW a (t, x,y). We use the notation 
P = (p x ,p y ) e R N+1 and W a (x,y, P) for the unique solution of (|2.5p . We have the following 
proposition: 

Proposition 2.2 (Barles, [3]). For any (t,x,y,P) € R x R 7 ^ 1 x R JV+1 , P = {p x ,p y ), the 
following estimates hold 

(i) m in (t,x, 2/ )eRxR' v + 1 ~F(t,x,y,P) < aW a {t,x,y,P) < m&x {t ^^ )mxR N+i -F(t,x,y,P); 

(ii) There exists a constant K\ > depending on \\F(t,x,y,p x ,Py)\\oo an d C*2 smc/i i/iai 

max W a - min W a < Jfi. 

Further properties of W a (x,y, P) are given in the following lemma: 

Lemma 2.3. For any (t, x,y,P) elx R^ 1 x R^ 1 the following estimates hold 

(i) a\DpW a (t, x, y, P)\ < C\, where C\ is introduced in (F2); 

(ii) \aW a (t,x,y, P) + F(P)\ < aK\, where K\ is introduced in Proposition\2J^ 

(iii) W a _(t,x,y,0) = 0; 

(iv) IpFHoo < C x . 

Proof. Let us fix Q £ R^ -1-1 . The Lipschitz continuity of F, i.e. (F2), implies that the function 
W(t, x, y) = W a (t, x,y,P + Q) satisfies 

W t + F(t, x,y,P + DW) + aW <d\Q\ 

and then, by comparison 

aW(t,x,y) < aW a (t,x,y,P) + C 1 \Q\. 
A similar argument shows that aW(t,x,y) > aW a (t,x,y,P) — C\\Q\. It then follows 

a\W a (t,x,y,P + Q)-W a (t,x,y,P)\ < d\Q\, 

which proves (i). 

Let us turn out to (ii). We claim that 

a := a max W a > -F(P). 

RxR N + 1 

Indeed, W a (t, x, y, P) is a supersolution of 

Wt + F(t, x,y,P + DW a ) = -ii. 

Let V be a bounded subsolution of (jl.7p . then by comparison between W a + fit and V — F(P)t, 
we have 

V(t, x, y) - W a (t, x, y) < V(0, x, y) - W Q (0, x, y) + t(F(P) + /x). 

Since V and W a are bounded, dividing by t > and letting t tend to +oo, we obtain fi > —F(P). 
Then from (ii) of Proposition EJ[ for (t, x, y) G R x R^ 1 , 

aFfi, x, y,P)>a min W Q > a max VT* - aKi > -F(P) - aK x . 
A similar argument shows that 

aW a (t,x,y,P) +F(P) < aKx; 

this concludes the proof of (ii). 

Property (iii) follows from F(t, x, y, 0,0) = and the uniqueness of the periodic solution of 

Finally, (iv) is an immediate consequence of 

F(P) - F(Q) < 2aKx + aWDpW^^P - Q\ 
and of (i). □ 
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We conclude this section by recalling some properties of the solutions vP and u e . 
Proposition 2.4. There exist constants Ct,L > such that for any (t,x), (s,y) G [0, T] x R^ 

(2.6) \u e (t,x)\, \u°(t,x)\ < C T , 

(2.7) \u (t,x)-u°(s,y)\<L(\t-s\ + \x-y\). 

Moreover, for any t G [0, T], the Lipschitz constant of u°(t,-) is the Lipschitz constant of the 
initial datum uq. 

Proof. By comparison 

\u e (t, x) — Uq(x)\ < C()t 

where Co = max.^ jd-ciwoh ^ \H(x,y,p)\. This implies (12. 6p for u e . Similarly can be showed the 
same estimate for u°. 

The Lipschitz continuity of u° follows from the comparison principle for (jl.2|) , see [2] , Theorem 
III.3.7 and Remark III. 3.8. □ 

2.3. Proof of the main result. This section is devoted to the proof of Theorem 12.11 We are 
going to show that for any T > 

sup \U e (t,x,y) - U°(t,x,y)\ < Ce T e^, 
[0,T]xR JV + 1 



where C does not depend on T. Since U e (t,x,y) = u e (t,x) — y and U°(t,x,y) = u°(t,x) — y, 
this estimate automatically gives f)2. 1 1) . 

Let us consider a function ^ : R ^ M with the following properties 



( i / / 



(2i 



s) > 0, for any s£l, 
lim (j)(s) = 1, lim (p(s) = 0, 

s — ^"l-OO s — — oo 

[ \<f>( 9 ) - X (s)\, W(s)\ < I g I> for any a G R, 



where we have denoted by x( s ) the heaviside function defined by 

X(s) = 

For n G N, e, 5 > 0, let us define the function 



1, for s > 0, 
0, for s < 0. 



^ 5 {s) := ^ e 0(^j-6(n+l). 

i=—n ' 

Then we have: 

Lemma 2.5. Assume (|2.8p . Then for any s G R ; the limit lim n ^ +00 ip™' S (s) exists and the 
function ip\: 

<p'(s):= Km ^ 5 (s) 
n— »+oo 

is o/ c/ass C 1 urai/i (yf)'(s) > for any s G R. Moreover 

(2.9) lim ^( S ) = { + ^f = ^. 

v ; 5 ^o+ \ *e, »/«e < s < (i + l)e. 

See the Appendix for the proof of the lemma. 
Let us define 

U e ' S (t,x,y) := V s € (W(t,x,y)). 
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Since F satisfies the "geometrical" assumption (F4), the function U e,s is solution of 

(2 1Q) f Uf + F (j , f , D x U e ' 5 , D y U^) = 0, (t, y) G (0, T) x R N +\ 

\ U e ' S (0,x,y)=<p 5 e (u (x)-y), (x,y) G R N+1 . 

By stability of viscosity solutions, see e.g. [7], the limit U e (t,x,y) of U e,s (t,x,y) as 6 — > + 
is a discontinuous viscosity solution of (|2.1U|) with initial datum <p £ (uo(x) — y), where (f t (s) = 
lim ( j_ > Q+ ip\ (s). This means that (U e )* = limsup^_ >0+ C/ e,<5 (resp. (t/ 6 )* = lim inf „ 5^0+ ^ e ' <5 ) is a 
viscosity subsolution (resp. supersolution) of (|2.10p . and (U e )*(0,x,y) < (ip t )* (uq(x) — y) (resp. 
(U e )*(0,x,y) > ((p e )*(u (x)-y)). Moreover, by §2M) 

{ie, if ie < U £ (t,x,y) < (i + l)e, 

(t - l)e + 0(O)e, if (t, x, y) G Int{£/ e = ie}. 

At the points (t,x,y) G d{U e = ie}, the value of U e depends on the lower semi-continuous or 
the upper semi-continuous envelope that we consider in the definition of discontinuous viscosity 
solution. In particular, since U e is continuous, U e has the following properties 

Wy{t,x,y)-U e {t,x,y)\, \(U%(t,x,y) - U e (t,x,y)\ < e for any (t,x,y) G [0,T] x R N+1 
and 

(2.12) DU e {t, x,y) = if U e (t, x, y) + ie, i G Z. 

Condition (|2.12|) implies that U € is actually a solution of 

j Ui + F(^,^,^,D x U%D y U £ ) =0, (t,x,y) G (0, T) x M Ar+1 , 
1 [7 e (0, x, y) = ^K(x) - y), (x, y) G M iV+1 . 

Indeed, when ie < U e (t,x,y) < (i + l)e, for some i 6 Z, the function C/ e is constant in a 
neighborhood of (t,x,y). Then the result follows from the fact that F(t,x,y,0) = 0. On the 

other hand, when U € (t,x,y) = ie, by periodicity, F (f , f , ^±a, p) = F(f,f,f,P). 

In order to estimate \U e — U°\ it is convenient to estimate |{7 e — ?7°|; indeed, — does not any 
longer appear in the equation satisfied by U e . 

Let us define V e (t,x,y) = e~ t U e (t, x, y) and V°(t,x,y) = e~ l U a {t,x,y). The functions V € 
and V° are respectively solutions of 

r91 o^ / V{ + V e + F{l2,\,D x V',D v V e )=0, (t,x,y)e(0,T)xR N +\ 

\ V*(0,x,y) = <p £ (u o (x)-y), (x,y) G R N +\ 

and 



(2.14) 



V t ° + V° + F(D X V°, D y V°) = 0, (t, x, y) G (0, T) x R^ +1 , 
F°(0,x,y) = Uo (x)-y, (x,y) eR N+1 . 

For alleviating the notations, let us denote a vector of R^" 4 ' 1 by X = (x,xn+i), where x G R N 
and xn+i G R. We first estimate from above the difference (V e )* — V°: for this, let us introduce 
the auxiliary function 

<&(*, X, S , Y) = (VT(t, X) - V°(s, Y) - eW a (-, -, 
(2 15) \ e e e 

\X-Y\ 2 \t-s\ 2 r 2 ?? 

2e^ 2^ l' ' T-t' 
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where a = e e , 9, (3, a, r, rj G (0, 1) will be fix later on and (3 and 9 satisfy 
(2.16) < 6* < 1 - /?. 

In view of $ZM), (|2TTT|h (i) of Proposition O and (|2~3l) . 



$(t,X, s,Y) <2C T + e + \x N+1 - y N+1 \ + -CV g— ^ - ' a - ~\X\ 



for all (t,X), (s,Y) € [0, T] x M iV+1 . Hence, <3? attains a global maximum at some point 
(7, X, s,F) G ([0, T] x M 7V + 1 ) 2 . Standard arguments show that 7, s < T for a small enough. 

Claim 1: There exists a constant Mi > independent of e such that < Mi(l + lyjv+il)- 
The inequality $(7,X,7,F) < $(7,X,s,F) and Proposition (JITIj) imply 

17 -si 2 



< v"(t,Y) -y u (s,y) < je-* - e - s ||c/ u (7,r)| + e - s |c/ u (t,y) -c/ u (s,y)| 



2(T 

< |i-s|(C T + lyjv+il) + 
from which Claim 1 follows. 

Claim 2: There exists a constant M2 > independent of e and T, such that -g 1 < Mi. 
The inequality $(7,X,s,X) < 3>(7,X,s,F) implies 

< V°(S,X) - V°fr Y) + eW« (I, f ,0) - eW a (I, f , ^- F% 
Using (|2.7[) , (i) of Lemma 12.31 and (|2.16p we then infer 

A ~ 5 ' < (L + 1)\X -Y\ + -Ci X ' = {L + l)\X -Y\+e 1 - e -< 3 C 1 \X -Y\ < M 2 \X -Y\ 



This concludes the proof of Claim 2. 

Claim 3: There exists a constant M3 > independent of e such that r|X| 2 < M3. 
The inequality $(7, 0,s, 0) < <3?(7,X,s,F) implies 

^|X| 2 < (F £ )*(t,X) - V°(a,F) + V°(s,0) - (V e )*(t,0) +eW a Q,0,o) - etT* 

Then, using (|23|) . (I2TTT) . Claims 1 and 2, (iii) of Lemma E3J (i) of Proposition and (1231) . 
we deduce 

r -\X\ 2 < e- J [U"(lX) - U\s,Y)) + \e- J - e~° \\U° (s ,Y)\ + e 
+ y°(5,0)-(^)*(7,0)-eH/ a f-,-, A * 



e e e 



< 4C T + M 2 e^ + |7 - s|(C T + \y N+1 \) + 2e + -Ci-^ 1 



a e 



< C + 2aM 1 \y N+1 \ 2 < C + 2aM l \X\ 2 , 



and Claim 3 follows by choosing a < 
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Now, suppose first that t = 0, then 



(Vr(t, X) - V°(t, X) - eW a I Oj - - "-\X\ 2 

< (^T(u (x) - x N+l ) - V°(s,Y) - eW a (o, 

for any (t,X) G [0,T] xf N+1 , from which, using (i) of Proposition E2J (hi) of Lemma E31 
and Claim 2, we deduce 

(Vy(t,X) - V°(t,X) < (^)*(u (x) - x N+1 ) - V°(s,Y) + + r -\X\ 2 + e 1 - e C 1 M 2 . 

Letting a, i] and r go to + and using (|2,lip and Claim 2 we obtain 
(Vr(t,X) -V°(t,X) < ( Ve y(u (x)-x N+1 ) - (u Q (y)-y N+1 ) + Ce l - d 

< (<p e )*(u (x) - x N+1 ) - (u (x) - x N+1 ) + {L + l)\X-Y\ + Ce l - e 

< C{/ + e 1 - 9 ) + e, 

which implies 

(2.17) U e (t, X) - U°(t, X) < Ce\cP + e 1 ' 6 ). 

The same estimate can be showed if s = 0. 
Next, let us consider the case t, s > 0. 

Claim 4: There exists a constant C > independent of e and T such that 

^ + 7^ + (VrilX) + F (1^1) < + e*). 

The function 

't X X-F\ |X-F| 2 r 2 |t-s| 2 7/ 



2(7 T-t 



(2.18) (t,X) - (V*)*fe *) " ^ f > " " ^1 

has a maximum at (t,X). By adding to <5 a smooth function vanishing with its first derivative 
at (i, X), we may assume the maximum is strict. 
Next, for j > 0, let us introduce the function 

^(t,s,X,Y,Z) : = (Vr(t,X) - eW a ( S ,F, - ~\\X\ 2 - 

!/ -(\t-es\ 2 + [X - Z\ 2 + |X -ell 2 ). 



T-t 2 

Let Pj = (tj, Sj,Xj,Yj, Zj) be a maximum point of on the set 

4 := B(t, l)xB 0, 1^ x B(X, 1) xB lj x B(X, 1). 

Since (I, X) is a strict maximum point of (|2.18p . tj —*-t, Sj — ► -, Xj, Zj — > X and lj — > — as 
j — » +oo. Then, for j large enough, Pj lies in the interior of A Moreover, standard arguments 
show that 

(2.19) j\tj-esj\ 2 , j|Xj-Zj| 2 , j|Xj-eljj 2 ^0 as j -»■ +oo. 
Remark that this implies in addition that 

(2.20) 2j\tj - eSj\\Xj - eYj\ < j\tj - esj\ 2 + j\Xj - eYj\ 2 ^ as j -> +oo. 

Since (y e )* and W^" are respectively viscosity subsolutions of f)2. 13j) and supersolution of (|2.5p . 
we obtain 
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tj — S 7] 



II ■ ■- 

(2.21) 



+ F [—, —, + rXj + j(Xj - Z 5 ) + j(Xj - eYj) ) < 



and 

(2.22) j( tj - es 3 ) + aW a ^Sj,Yj, + F (sj, Yj, + j(Xj - eYjfj > 0. 

Subtracting (|2.2ip and (|2.22p and using the Lipschitz continuity of F, assumption (F2), we get 
(2.23) 

tj p + (^)2 + <rnti,xj) - ^ (s jt Y jt ^=pj 

< Y (\tj - e Sj \ + \X 3 - eYj\) ( |ZJ j y| +j\Xj - eY^ + C x ~ Zjl + r\Xj\ + j\X 3 - Z 3 

Let us estimate j \Xj — Zj\. From the inequality tyj (tj , Sj , Xj ,Yj,Xj) < (tj , Sj , Xj , Yj , Zj ) we 
deduce that 

J -\Xj - Z,| 2 < eW° fa, Yj, " W a (sj, Yj, 

and using (i) of Lemma 12.31 we get 

i\Xj - Zj? < c x ^i^ = c^-«-e\Xj - Zj\. 

Then 

(2.24) j\Xj -Zj\< 2C 1 e l - e ~ 13 . 

Then, passing to the limsup as j — > +oo in (|2.23[) and taking into account Claim 2, (|2.19p and 
(j23Uj) . we obtain 

(2.25) lp + + (Vy(t,X) - aW« (?, * ^) < C(e^ + r\X\). 



1 



By Claim 3, r\X\ < r^M^, hence choosing r > such that r^M^ < e l ~ e ~@, we have r\X\ < 
f i " \ 

Finally, Claim 4 easily follows from (|2.25p . Claim 2 and the following inequality 



which comes from (ii) of Lemma 12.31 . 

Claim 5: There exists a constant C > independent of e and T such that 



The function 

'£ X X-y\ |X-y| 2 |i- s | 



(a, Y) - 0( S , y) := F°( S , y) + eW* f , ^) + + 



2(7 
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has a minimum at (s, Y), consequently (0, 0) € D </>(s, Y). If we set 

V(s,Y) := V°(s,Y) + \^pL + ILiE, W{Y) := eW<* (I, * A " V 



by properties of semijets of Lipschitz functions, see e.g. Lemma 2.4 in [41], there exists Q G f Ar+1 
such that 

(0, Q) G ZrVftF) = D-V°(s,Y) - (^-, ~ Q G D~W{Y). 

Since y° is a supersolution of (|2,14p . we have 

(2.26) L^ + yO^+^^I + Q^ > . 



(7 



By (i) of Lemma 12. 3\ 



a e p 



from which we get the following estimate of Q: 

(2.27) \Q\ < de 1 - -! 3 . 

Then, Claim 5 follows from (|2.26p using estimate ()2.27|) and the Lipschitz continuity of F as- 
sured by (iv) of Lemma [ 



Claims 4 and 5 imply 

(V e )*(t,X) - V°(s,Y) < C{e l -°- p + e 9 ), 

for some constant C independent of e and T. Since (t,X,s,Y) is a maximum point of <3?, we 
have 

(V e )*(t,X) -V°(t,X) < <S>(t,X,s,Y) + eW a (-,—,o) + ^\X\ 2 + '' 



e e 



T-t J 



for all (t,X) G [0, T] x R N+1 . Then, by (iii) of Lemma 

(vy(t,x) - v°(t,x) < (vy(t,x) - v°(s,Y) - e\v a Q, |, *-zTj + r .\ X \* + ^ 

< C(e^-P + e°) + -Cj*^ + r \X\* + -^- 
a 2 1 — t 

<C{e^ + e°) + r -\X? + ^- t , 

for some positive constant C. Hence, sending r, r/, — > + and taking into account (|2.1ip . we get 

U"{t, X) - U°(t, X) < Ce\e x - Q - (i + e e ). 

Then, from the previous estimate and (12. lTf) . we can conclude that for all /3, 9 G (0, 1) satisfying 
(I2.16P we have 

U e {t, X) - V°(t, X) < Ce'ie 1 - - 13 + e e + e?), 
for all (t, X) G [0, T] x M Ar+1 . The optimal choice of the parameters is 9 = f3 = |, which gives 

sup {U e (t, X) - U°(t, X)) < C(h . 

[0,T]xR N + 1 

The opposite inequality follows by similar arguments, replacing (V 6 )* with V° and V° with 
(V e )* in (|2.15p . and the proof of Theorem 12.11 in the general case is complete. 
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Now, let us consider the case when uq is affine. Let us suppose that uq(x) = p ■ x + cq for 
some p G 1^ and c G K. In this case, the solution of (|1 .2|) is u°(t,x) = p ■ x + cq — H(p)t. Let 
V be a bounded viscosity supersolution of (II. 7p with p x = p and p y = —1. Let us define 



V e (t,X) = U°(t,X)+eV 



t X 



oo • 



Since «o(x) — y > ip e (uo(x) — y) — e then V e (0,X) > ip e (uo{x) — y) — (M + l)e where M = ||V| 
Hence, it is easy to check that V e is a supersolution of 

Vi + F (I, ^,D x V e ) = 0, (t, X) G (0, T) x R Ar+1 , 

V*(0, X) = 99 e (n (x) —y) — (M + l)e, (x, y) G M Ar+1 . 

By comparison we get V e (t,X) > (U e )*(t,X) - (M + l)e and this implies that U°(t,X) - 
U e (t,X) > —Ce. A similar argument shows that U°(t,X) — U € (t,X) < Ce and this concludes 
the proof of the theorem. □ 



3. Approximation of the effective Hamiltonian by Eulepjan schemes 

In this section we give an approximation of the effective Hamiltonian F(P). To this end, we 
introduce an approximation scheme for the equation (|2.5p and for simplicity we only discuss the 
case N = 2. Given Nx and Nt positive integers, we introduce At = 1/Nt, h = 1/Nx and 

R l : = = ( x i,Vj) I x i = ih, yj = jh, i,j G Z}, 

^■At ■= {t n = nAt\n G Z}. 
An anisotropic mesh with steps hi and hi is possible too; we take hi = hi only for simplicity. 
We denote by W i '• ' our numerical approximation of W^ 3 ' 01 at {tnt Xi? Uj 

) G M A t x R|. For (1231) 

we consider the implicit Eulerian scheme of the form 



w n+l,P,a _ W n,P,a 

(3.1) -a At " J +aW^ 1 ' P ' a + S(t n ,x i ,y j MW n+1 ' P ^, j ) = 0, 

where 

S(t n ,Xi,yj,h, [W]u) 

(3 2) 

= K^iZi 9 yi>( A ^)ij+P*,(AfW)i-ij+Px 



and 

We make the following assumptions on g: 



Wi+i j ~ W id _ Wjj+i - w itj 

h ' w /hi ~ I 



(gl) Monotonicity: g is nonincreasing with respect to its fourth and sixth arguments, and 

nondecreasing with respect to its fifth and seventh arguments; 
(g2) Consistency: for any t G R, (x,y) G R 2 and {q x ,q y ) G R 2 

g(t,x,y,q x ,q x ,q y ,q y ) = F(t,x,y,q x ,q y ). 

(g3) Periodicity: for any t G R, (x, y) G M 2 and Qel 4 

#(i + 1, x + 1, y + 1, Q) = g(t, x, y, Q); 

(g4) Regularity: g is locally Lipschitz continuous and there exists C\ > such that for any 
t G M, (x,y) G M 2 and Q G M 4 

|£>0^(t,a;,i/,Q)| <C i; 
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(g5) Coercivity: there exist C 2 , C3 > such that for any i G R, (x,y) G R 2 , (91,92) G R 2 



1 



g(t,x,y, 92, 0,0) > C 2 (|^| 2 + |?2"| 2 )* - ^35 

(g6) For any t G R, (x, yi), (x, y 2 ) G R 2 , 91, 92 G R 

g{t,x,y 1 ,q 1 ,q 2 , 0,0) = g(t, x, y 2 , 91, 92, 0, 0). 

The points (gl)-(g4) are standard assumptions in the study of numerical schemes for Hamilton- 
Jacobi equations. The coercivity hypothesis (g5) can be substituted by the weaker condition 

lim g(x, y, q 1 ,q 2 , 93, 94) = +00 
<?i +<? 2 -^+°° 

if g (and hence F) does not depend on time. If g is homogeneous of degree 1 w.r.t. Q, then the 
two coercivity conditions are equivalent. 

As an example, we suppose that the Hamiltonian F is of the form F(t, x, y,p x ,P y ) = a(t, x)\p x \+ 
b(t,x,y)\p y \, with a and b Lipschitz continuous functions and a(t,x) > C 2 > 0; we consider a 
generalization of the Godunov scheme proposed in |15| : 

g{t,x,y,q 1 ,q 2 ,q 3 ,q A ) 

= a(t,x)[(9f) 2 + (9 2 + )¥ +b + (t,x,y)[(q^ f + (g 4 + ) 2 ]^ - b~ (i, x,y)[(qj) 2 + (q^ . 

where q + = max(9, 0) and q~ = (—q) + - Then hypothesis (gl)-(g6) are satisfied. 

The following theorem is the discrete version of the analogous result in [3] for the exact 
solution W p ' a of (1231). 



Theorem 3.1. Assume (gl)-(g6). Then we have 

(i) For any P = (p x ,p y ) G R 2 , a, h, At > there exists a unique {W^j P ' a ) periodic solution 

(ii) There exists a constant K\ depending on \\F(-, •, ■, -P)||oo; Ci ^ n (g4)> C 2 , C3 in (g5), p x 
and p y , but independent of a, h and At such that 

max W ij ■ - mm ^ < K\\ 

i,j,n " 

(iii) There exists a constant F^ t (P) such that 

(3.3) YimaW?f> a = -Fi\P) Vi,j,n; 

a— >0+ 

(iv) F^ t (P) is the unique number A^* G R such that the equation 

W n+1 ' p - W n ' P _ A , 
(3-4) l ' J — M + S(t n , x hyj , h, [W n+l > p \ id ) = Af 

admits a bounded solution. 

Proof. A proof of the existence of a unique solution of (|3,ip in the uniform grid on the torus 
with step h is given in [6]. 

Let us prove (ii). First, remark that by comparison with constants we have 



(3.5) m&x\aW"f> a \ < Q 



0, 



where Cq := \\F(-, •, •, -P)||oo- Next, let us define 



TTT 71 TTTll,P,a 

W \ := max W i ] ' . 



We claim that W4 satisfies 

1 yFTTt 

' + aW-^ + S(t n , Xi , h, [W'^i) < 0, 



777™+ 1 777™ 

W i ~ W i . — n+l , —,, , r — n+1-, 



At 
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where 

S(t n ,Xi,h, [W]i) := mm g(t n ,Xi,yj, (AfW)i + p x , (AfW)i-i + p x ,p y ,p y ). 
j 

Indeed, for any i and n, denote by jn in \ the index j such that W i = max,,- 
then 

w n+l,P,a _ yyn^a w n+l,P,a _ w n,P,a +1 

l >i(i,n+l) M(i,n + 1) > M(i,n+1) l >i(i,n) _ W j _ ^ i 



At ~ At At 

p^n+l,P,a _ y^n+l,P,a y^n+l,P,a _ yyn+l,P,a 
I A +TT/n+l,P,a\ _ _ t+1 'J(i,n+l) M^n+l) , H-l,jj + i, n+ i ^%n+l) _ / A +t77» 1 + 1 n, 

y^n+l,P,a _ ^yn+l,P,a ^n+l,P,a _ yyn+l,P,a 

(A + W n+1 ' P ' a ) - — ' ,:)(i '° +1> »- 1 J(i,n + l) > M(i,n + 1) »-lJ(i-l,n+l) _ (A^+pp^+V i 

and 

yyn+l,P,a _ ^n+l,P,a 

Since (W,?J ,a ) satisfies (IBTTj) . using the monotonicity assumption (gl), we get 
^ + A ~^ + aWT 1 + 5(t n , ft, W n+ \) 

At l>3(i,n+l) 

+ g{t n ,x i ,y- j , (A^" +1 )i + p x , (A^W n+1 ) i - 1 +Px,Py,Py) 
w n_+l,P,a _ jyn.P.a 

~ At 



(A ^n + l,P, Q)i _^ +i) + ^ (AtW^^^+Py) 



<o, 

as desired. Then, by (g4), (g5) and (|3.5|) . we see that W™ satisfies 



+ 1 777™ . 1 



+ C 2 \[(AJW )i +p x }~\ + {[(A+W'^)^ + Px ] + \ 2 - < 



At 

where K\ = Co + C3 + 2C\ \p y \ . In particular we infer that 

W^ +1 -Wl< K x At, 

which implies that if n > m then 

(3.6) -W™ <K 1 (n-m)At = K 1 (t n -t m ). 

Next, let us consider 

Wi = max W™. 

n 
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Similar arguments as before show that W \ satisfies 



C 2 (\[(A+W)i+ Px ]-\ A + |[(A+W)i_i + Px ] + \ 2 ) < Kx, 

which implies the existence of a constant K 2 > depending on Co, C±, C 2 , C%,p x and p y such 
that 

(3.7) max|(A+W)i| < K 2 . 

i 

Now, let (ii,ni) and (i 2 ,ra 2 ) be such that maxj^W™ = W? 1 and mhij n = W™ 2 , and let 

rij 2 be such that PFj 2 = max n VF™ = W™* 2 . By periodicity, we may take — :rj 2 | < 1 and 
< t Hi — t n2 < 1. Then using (|3.7p and (|3.6p . we get 



< Wi 2 + ^l^ii - x i2 \ 

rn 2 



<W l2 2 +K 1 (t ni -t n2 ) + K 2 



12 

r«2 



''2 

Then we have proved that 

(3.8) max - min VF™ < if , 

where Kq depends only on Co, Ci, C2, C3,p x and 

Next, we consider the behavior of Wj^ ,a in j. We claim that 

^ijf '° + p y y n ^ ^Cjf ' a + p y y v if - J' 2 and ^ < °' 

(3.9) W^f' a = W4 P ' a for any * , j 2 if p y = 0, 

Wj5f + ?Wi > Wjjf '° + W2 if ii > h and pj, > 0. 
Let us consider the case p y < 0. Suppose by contradiction that 

M := . max. (W?£ a - W^ a + p y (y n - y n )) = Wlf> a - W^+p^ - y ) > 0. 
Then j l > j 2 + 1. We have the following estimate 

W n,P,a _ W n,P,a W %P,a _ w %P,a 

< 0. 



W n ' p ' a 




i,h 


W-> P '- 




W-i P ' a 




h 






h 








W n,P,a 


w %p, a 




W n z P,a 




h 






h 





Similarly 
and 



Mi hh > M2 

At ~ At 
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Moreover, we have 

rn,P,a 



W n,P,a 




W--' a 




h 








w?± p > a 


h 



similarly 

(A+^"' p ' a ) ?Ji _ 1 + Py > 0, (A+W^)^ > 0, (A+W*> p < a y i32 _ 1 + Py < 0. 
Then, since W?f> a satisfies ([33]), using assumptions (gl) and (g6), we get 

+ g(t m xj, y- ]2 , (AtW*> p '«)- i3i + Px , (A+ ^ + p x , 0, 0) = 0. 

This implies that 

< «M = a(W^ a - + Py ( yji - yjj) < ap y {y- h - ^)< 0, 

which is a contradiction and this concludes the proof of (|3.9p for < 0. The case p y > can 
be treated in an analogous way. 

Now, to prove (ii), we use the properties (|3.8p and (|3.9|) of W™' P,a and again we only consider 

rni,P,a _ w n,P,a 



the case p a < 0. Let («'i,ji,ni) and {11,32^2) be such that W^^' 01 = maxj jJ>n ,a and 

»2,J 



W^ 2 j 2 = minij i7 i W™j P,a . Let j be such that W™ 2 = W^ 2 - P,a . By periodicity, we can take 



< Uj — yj 2 < I and — a?j 2 | < 1. Then 



y^ni ,P,a ^1 



< ^r 2 2 + ^0 

<^f' a +^-yj)+^o 

<^f' a -p.+^o, 

and this concludes the proof of (ii). 

The property (iii) easily follows from (ii) and (|3.5p . Indeed, from (|3.5|) . up to subsequence, 

a mirLij^ W™' P,a converges to a constant — F h (P) as a — > + . Then from (ii), for any i,j,n, 
we get 

|aWij +F h (P) | < \ammW i ' +F h {P)\ + a\WA -mmW- \ 

' J Mi n i>,3,n 

< \a min W?f' a + F^\ + aK x -► as a -» + , 

and (iii) is proved. 

Let us turn to (iv). Let us define Z™j ' a = W^j ' a — minj Jin W™j ,a . By (ii), up to subsequence, 
(Z™' P,a ) converges to a grid function (Z^ ,P ) as a — > + . The grid function (Z^' P,a ) satisfies 

™n-j-l,P,a ^n,P,a 

l ' 3 ~ l ' J + nz;;; + S(t n , Xi , yj , h, [Z n+1 > P >% tj ) = awinW?f' a . 

Letting a — > + , since by (ii) ,a ) is bounded and o: minj J ' ► F^ , we see that 

(Z™ ,P ) is a solution of (|3.4p with A^* = 
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To prove the uniqueness of a solution (A^*, {W™' P )) of (|3.4p . we show that if there exists a 

subsolution {Uj 1 ^) of (|3.4j) with A^* = Ai and a supersolution (V^' P ) of (|3.4j) with A^* = A2, 
then A2 < Ai. 

Let M = max,,,.!,/// - V?f) = V%£ - V%£. Then 

Tj-no,P _ rj-no — ljP y n o>P _ y"o— l,f 

Ai - At 

(A+c/^) lWo < (A+y^. p ) i0)io; (A+c/"o^) io _ 1)io > (A+y«'^) JO _ 1Jo; 

(A+C/"°' p ) lWo < (A 2 + F"°< p ) Wo , (AtU n °> p ) i0>j0 ^ > (A^V n °' p )i 0> j -i. 
From the monotonicity of g, 

jjTl ,P _ jjU -1,P 

Ai> i0 ' j0 At ioJo +ff(i no ,^ , %0 ,(A+^°Oi Jo+^,(A+^°- p ) io - 1 , io +^, 



(A+C7"»- p ) l0J0 +p„ (A+L^.- p ) io>J - _ 1 



yn ,P _ yU -l,P 

(A+y^> p ) JOiio + P „ (A+y^ p )^o-i 

> A 2 . 

This concludes the proof of (iv) . □ 

We need a more precise estimate on the rate of convergence of aW^^ ,P to F^ t (P): 
Proposition 3.2. Assume (gl)-(g6). Then for any i, j,n 

\aW^ P + Ft\P)\<K ia , 
where K\ = K\(P) is the constant in (ii) of Theorem \3.1[ 

Proof. As in the proof of (ii) of Lemma 12.3} the result follows from the comparison principle 
for ([53]) and (ii) of Theorem EQ □ 

Now, we are ready to show that the function is actually an approximation of the effective 
Hamiltonian F. 

Proposition 3.3. Assume (gl)-(g6). Let F^ be defined by (|3.3p and let F be the effective 
Hamiltonian. Then, for any PgM 2 

lim F$\P)=F(P) 
(At,/i)-»(0,0) 



uniformly on compact sets of 



p2 



Proof. To show the result we estimate W P ' a (t n ,Xi,yj) — W™j ' a . To this end, following the 
same proof as in [8] and pQ, we assume that 

sup \aW P ' a (t n , Xi , yj ) - aW^ P ' a \ = sup(aW P ' a (t n , x h Vj ) - aW?f' a ) = m > 0. 

i,j,n i,j,n 

The case when sup^ \aW p > a (t n , Xi, yj) - aW^ P ' a \ = sup i j >n (aWl i j P ' a - aW p > a (t n ,Xi, yj )) is 
handled in a similar manner. 

For simplicity of notations we omit the index P. Let us denote Wff At (t n , Xjj) := Wj 1 ^ 01 , 
(t n ,Xij) G R At x R|. For (X,Y) € R 2 x R 2 h and (t,s) £Rx R At , consider the function 



* ( / . A\ s . V ) - - n 1 V ' ( / . X ) - oil 7/. A/ ( > V ) + ( 5C + ^) (3 e (t - S) X - Y), 



18 



YVES ACHDOU AND STEFANIA PATRIZI 



where, as before, Cq = \\F(-, •, •, P)||oo and (5 t = f3 (|, ^) with (3 a non-negative smooth function 
such that 



f3(t, X) = 1 - \X\* - \t\ A , if \X\ Z + |t| 2 < i, 



if | < |A~| 2 + |t[ z < 1, 
if |A1 2 + Itl 2 > 1. 



/3 = 0, 

We have the following lemma: 

Lemma 3.4. The function ^ attains its maximum at a point (to, Xq, so, Yq) such that 

1' 

5 ' 

For the proof, see Lemma 4.1 in [8]. 
Lemma 13.41 (ii) implies that 

(3 e {t - s ,X -Y ) = 1 
Then, from the inequality ^(sq, Yq, sq, Yq) < ^(to, Xq, sq, Yq) we deduce that 



(i) ^(t ,X ,s ,Y ) > 5C + |m; 

(ii) p e (t -s ,X -Y ) > 3 



X -Y 


2 


to - s 


e 




€ 



(3.10) 



5C + y 



+ 



*o - s 



< aW a (t ,X ) - aW a (s ,Y ) < 2C . 



This implies that \t — sq\ — >• and \X — Yq\ — > as e — > 0. Moreover, since iy a and W^ At 
are periodic, we can assume that (to, -Ao, so, Yq) lies in a compact set of (R x M 2 ) 2 . Hence, from 
(|3.10p and the continuity of W a we get that 



(3.11) 



Xo-Y 



+ 



to - so 



as e — > 0. 



Since (to, A~o) is a maximum point of (t, X) — > aiy Q (t, X) + (5Co + y) /? e (t — so, AT — Yq), we 
have 



5C + f 



(3.12) 



o 



dtPefa -s ,X - Y ) + aW a (t , X ) 



+ / j /„. A'o, - 5Co + m/2 ^x/?.(to - so, X - Y ) + P ) < 0. 



n 



Let iojio and no be such that Xi j = Yq and so = t no . Since (so>^o) is a minimum point of 
(s, Y") a^" A( (s, Y") - (5C + m/2) &(t -s,X -Y), we obtain 



]4/ n 0,a _ w n o,a > 5Cq + m/2 



io+l>jo »oj'o 



o 



[/3 e (t - so, A~ - Y" - /iei) - /3 e (to - s , X - Y )], 



where e\ = (1,0) T . From the monotonicity of g, 



(3.13) At + aW ^o +9 I so, y 



5C + m/2 



a 



(A|/3 e (t - s ,X - -))io,io 



(A+W"°'«) io _ 1)jo + p x , (A+W n °' a )ioJo +Py, (A$W n °> a )io,Jo-l +Py) > 0. 



But 



h, 



\(A+(3 e (t - s ,X - -)) iodo - ei • Dy/3 £ (t - s ,X - y )| = ^|e^ Y /3 e (t - s ,X - y) ei |, 

for some Y belonging to the segment (Yq, Yq + hex). Assuming h small enough, so that Lemma 
E3(ii) implies that |t -s | 2 + |Ao-yo| 2 + /i 2 < 4> we obtain that D%r Y (3 e (t -SQ,X -Y) = Jr/, 
then 



(3.14) 



\(Af/3 e (t - s , A'o - -))<oJo ~ e i ' D rPe(to - so, A'o - io)] = 4- 
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Now, (|3.13p . (|3.14p and the monotonicity of g yield 

r no— l,a 



w no,a _ w <<-o- 

At 



+ aW%£ + g[s ,Y 



5C + m/2 



o 



ei • D Y /3 e (t - s ,X - Y Q ) +pz 



(AtW n °' a )i -ljo +Px, (AtW n °n*o,io +Py, (AtW n ^ ,j -i+Py) + Cih 



5C + m/2 



e 2 a 



> 0. 



Repeating similar estimates for the other arguments in g and for the derivative with respect to 
time, we finally find that 

5C + m/2 . 



(3.15) 



a 



-d a p e (to - so, Xq - Y ) + aW? Q °j«+ 



F(s ,Y 



5C + m/2 



a 



D Y (3 e (t -s ,X -Y )+P) +C 



,h + At 

e 2 a 



where C is independent of h, At, e and a. 

Subtracting (|3.12p and (I3.15|) and using (F2) we get 



(3.16) 



aW a (t ,X ) - aW? At (s ,Y ) < C 



h + At C 



e 2 a 



+ 



o 



X -Y 



+ 



C 



> o. 



*o - so 



n 



where C is independent of h, At, e and a. 

Choose e = e(At,h) such that e -» as (At, /i) -> (0,0) and -» as (At, h) -> (0,0). 

From (i) of Lemma 13.41 

sup |aW^>(t n , si, % ) - a^ P ' a | = m < sup * - (5C + y) &(t - s , X - Y ) 



i,3,n 



Then from (|3.16p and (|3.1ip . we obtain 



aW a (t ,X ) - aW? At (s ,Y 



C 



sup \aW^ a (t n , Xi , yj ) - aW?f' a \ < -o(l) as (At, h) (0, 0) 



From the previous estimate, (ii) of Lemma 12.31 and Proposition 13.21 we finally obtain 



\F(P)-F*\P)\ <K 1 a + K 1 a + -o(l), 

a 



and letting (h, At) — > (0, 0), we find that 



limsup |F(P) - F h (P)| < iTia + ifia, 
(At.h)— (0,0) 



for any fixed a > 0. This implies that hm(At,k)-i-(o,o) ^"ft = F(P)- Since = K\(P) and 
= K\(P) are bounded for P lying on compact subsets of M 2 , the convergence is uniform on 
compact sets. 



□ 



Remark 3.5. If F is coercive, then we can get an estimate of the rate of convergence of F h 
to F. Indeed, we have: 



At 



=rAi 



Fh - F \<(h + At) 



see Proposition A. 3 in pQ. 



At 



We conclude this subsection by recalling the principal properties of F h . 

Proposition 3.6. Assume (gl)-(g6), (Hl)-(H^). Then the approximate effective Hamiltonian 
F h is Lipschitz continuous with a Lipschitz constant independent of h and At and for any 

Fi\p x ,0)>C 2 \ Px \. 
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Proof. For the proof of the Lipschitz continuity of F, see the proof of Proposition A. 2 in pQ. 

Let us show the coercivity property. Let (W™j be a solution of (|3,4|) for P = (p x ,0). Let 
(zo>jo ; n o) be a maximum point of (W™f' a ), then 

^n 0; P,a _ yyn -l,P,a 

zoJo At '°' J ° > 0, (A+W n °' P ' a ) io , jo < 0, (A+W n °> p >«) io - 1 j > 0, 

{A + W no,P n _ Q < Q) {A + w n ,P n ^ , o _ i > _ 

By the monotonicity assumption (gl) and (|2.4p . we have 

F h\Px,0) > g(t no ,x io ,y io ,p x ,p x , 0,0) = F(t no ,x io , y io ,p x ,0) > C 2 \p x \. 

□ 

3.1. Long time approximation. A different way to approximate the effective Hamiltonian is 
given by the evolutive Hamilton-Jacobi equation 

V t + F(t,x,y,p x + D x V,p y + D y V) = 0, (t,x,y) G (0, +oo) xR N +\ 
V(0,x,y) = V o (x,y), (x,y) eR N +\ 

where Vb is bounded and uniformly continuous on M Ar+1 . Indeed, it is proved in [3] that (|3.17j) 
admits a unique solution V which is bounded and uniformly continuous on [0, T] x M Ar+1 for any 
T > 0, and satisfies 

lim V( t> X > y) =-F(P). 

t^+co t 

We approximate (|3.17p by the implicit Eulerian scheme 

yn + l,P _yn,P 



(3.17) 



(3 18) '' J At l,j + S(t n ,x i ,y j ,h,[V n+1 > p ] i , j ) = 

V i, =Vo(x i ,y j ), 



where S is defined as in (|3.2p . A proof of the existence of a solution V = (V^™' P ) of (|3.18p is 
given in [6] under assumptions (gl)-(g5). 

Let W = (W™J P,a ) be a solution of (|3.4p . then by comparison, there exist constants c and c 
such that 

£ + Kf' a ~ nF%\p)At < Vlf < c + W?f> a - nFf (P)At. 
Since VF is bounded, this proves that 

V n -: P a x 

lim -M_ = _Ff(P). 

3.2. Approximation of the homogenized problem. We now come back to the A^-dimensional 
homogenized problem (II. 2p . From Theorem 11.31 we know that if H is the effective Hamiltonian 
in (jl.2p . then -ff(p) = F(p,—1) for any p G M^. Hence, from Proposition 13.31 the discrete 
Hamiltonian 

H£\p):=F£\ P ,-l), 

is an approximation of H(p) for any p G R^. 
As in pQ, we approximate (jl.2p by the problem 



(3.19) 



^«At,h + H h (Du Atth ) = 0, (t, x) G (0, +oo) x 
UAt,h(Qi x ) = uq(x), x G R N , 



where h and At are fixed, and u$ is the same initial datum as in (|1.2I) , 

By Proposition 13.61 H h is Lipschitz continuous and coercive, so (|3.19p has a unique viscosity 
solution UAt h which is an approximation of the solution u° of (jl.2p : 
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Proposition 3.7. Let u° and «At,fi, be respectively the viscosity solutions of (|1.2|) and (|3.19p . 

T/ien /or any T > 

(3.20) sup \u At ,h -u°\^0 as (At, h) -» (0, 0). 

[O.TJxR^ 

Proof. If Lo is the Lipschitz constant of the initial datum uq, then, by Proposition 12.41 the 
functions u (t, •) and u A t,h(t,-) are Lipschitz continuous with same Lipschitz constant Lq. By 

Proposition 13.31 the approximate Hamiltonian converges to H uniformly for \p\ < L. Hence 
(|3.20p follows by the following proposition, which is a standard estimate in the regular pertur- 
bation theory of Hamilton- Jacobi equations (see Theorem VI. 22.1 in [2]) 

Proposition 3.8. // there exists ij > such that if Hi, i = 1,2, satisfy (H1)-(H3) with 

ll-Hi - H 2 \\oo < i], 

and if Ui, i = 1,2, are viscosity solutions of 

u t + Hi{Du) = 0, (t, x) G (0, T) x R N 
u(0,x) = u (x), x€R N , 

where uq is bounded and uniformly continuous on R N , then, for some constant C, 

\\ui - u 2 \\oo < Crj. 

□ 

Remark 3.9. In order to compute numerically the approximation of vP, we need further dis- 
cretizations. Indeed, we have approximated H(p) by H h (p) for any fixed p € R . Since it is 
not possible to compute H h (p) for any p, one possibility is to introduce a triangulation of a 
bounded region of and compute H h (pi), where pi are the vertices of the simplices and to 
approximate all the other values H h (p) by H hk (p), where H h k is the linear interpolation of 



H h and we denote by k the maximal diameter of the simplices. The solution u^ t h of 
(3.21) 



9tu k Atih + H^ k (Du k Atth ) = 0, (t, x) e (0, +oo) x R N , 
u k Ath (0,x) = u {x), xeR N , 



is an approximation of u A t,h as k — > and hence, by Proposition 13.71 °f u ° as (At, h, k) — > 
(0,0,0). Finally, discretizing (|3,2ip by means a monotone, consistent and stable approximation 
scheme, we can compute numerically an approximation of the solution u° of 11.21 See [1] for 
details. 

4. Numerical Tests 

The present paragraph is devoted to the description of numerical approximations of the ef- 
fective Hamiltonian. 

4.1. Results. 

4.1.1. First case. We discuss a one dimensional case where the Hamiltonian is 

H(x, u,p) = 2cos(2-7rx) + sin(87ru) + (1 — cos(67rx)/2)|p|. 
We have used two approaches for computing the effective Hamiltonian. 
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(gl) Barles cell problem: the first approach consists of increasing the dimension and consid- 
ering the long time behavior of the continuous viscosity solution w of 

, , w t + F(x,y,p + D x w,-l + D y w) = 0, (t, x, y) £ (0, oo) x R x R, 

1 ' ' w(0,x,y) = 0, (x,y) G R x R, 

where F is given by (|1.4|) . In the present case, from the periodicity of H with respect to 
x and u, w is 1-periodic with respect to x and 1/4-periodic with respect to y. We know 
that when t — > oo, w(t, •, -)/t tends to a real number A and that H(p) = —A. 
For approximating (|4,ip on a uniform grid, we have used an explicit Euler time marching 
method with a Godunov monotone scheme (see [SI [16] ) . A semi-implicit time marching 
scheme which allows for large time steps may be used as well, see PQ, but very large time 
steps cannot be taken because of the periodic in time asymptotic behaviour of w. 
Alternatively, we have also used the higher order method described in [13], see also |14j . 
It is a third order TVD explicit Runge-Kutta time marching method with a weighted 
ENO scheme in the spatial variables. This weighted ENO scheme is constructed upon 
and has the same stencil nodes as the third order ENO scheme but can be as as high as 
fifth order accurate in the smooth part of the solution. 
(g2) Imbert-Monneau cell problem: when p is a rational number {p = ^), instead of consid- 
ering a problem posed in two space dimensions, one possible way of approximating the 
effective Hamiltonian H(p) is to consider the cell problem 

, . „s v t + H(x,v + p ■ x,p + Dv) = 0, (t, x) € (0, oo) x R, 

( " ' v(0,x) =0 i£R. 

This problem has a unique continuous solution which is periodic of period q with respect 
to x (in fact, the smallest period of v may be a divisor of q). From [11] (Theorem 1), we 
know that there exists a unique real number A such that V ( T,X > converges to A as r — > oo 
uniformly in x, and that H (p) = —A. Moreover, when t is large, the function v(t, x) — Xt 
becomes close to a periodic function of time. In what follows, (|4.2p will be referred to as 
Imbert-Monneau cell problem. Note that the size of the period varies with p and may 
be arbitrary large. This is clearly a drawback of this approach which is yet the fastest 
one for one dimensional problems and moderate values of q. 

For approximating (|4.2p on a uniform grid, we have used either the abovementioned 
explicit Euler time marching method with a Godunov monotone scheme or the third 
order TVD explicit Runge-Kutta time marching method with a weighted ENO scheme 
in the spatial variable. 

In Figure [H we plot the graph of the effective Hamiltonian computed with the high order methods 
and both Imbert-Monneau and Barles cell problems. For Barles cell problems, the grid of the 
square [0,1] x [0,1/4] has 400 x 100 nodes and the time step is 1/1000. For Imbert-Monneau 
cell problems, the grids in the x variable are uniform with a step of 1/400 and the time step is 
1/1000. The two graphs are undistinguishable. It can be seen that the effective Hamiltonian is 
symmetric with respect to p and constant for small values of p, i.e. \p\ < 1.3. The points where 
we have computed the effective Hamiltonian are concentrated near 1.3 where the slope of the 
graph changes. Our computations clearly indicate that the effective Hamiltonian is piecewise 
linear. 

In order to show the convergence of V ( T,X ) and W ( T > X ) ; we take p = 1.3 so the space period of 

the Imbert-Monneau cell problem is 5. In Figure [21 we plot (left) and '^ r " (right) as 

a function of r, where (v(r)) is the median value of v(r, ■) on a spatial period. Both functions 
converge to constants when r — > oo and the limit are close to each other (the error between 
the two scaled median values is smaller than 10 -3 at r ~ 60 and we did not consider much 
longer times). In Figure [31 we plot the graphs of the functions w(t, 0, 0) — (w(t)) (left) and 
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4.5 



3.5 



2.5 



" Imbert-Monneau" — r 1 — 
"Barles" 




Figure 1. First case: the effective Hamiltonian as a function of p obtained with 
both Barles and Imbert-Monneau cell problems. 
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Figure 2. First case, p = 1.3. Left: the median value of w(t, -)/t on a period 
as a function of r. Right: the median value of v(t,-)/t on a period as a function 
of r 



v(r, 0) — (v(t)) (right). We see that these functions become close to time-periodic. In Figured] 
(top), we plot the contour lines of the function w(t, x, y)/r as a function of (x, y) for r = 60. In 
the bottom part of the figure we plot the graph of y — > w(t, 0.13, y)/r for the same value of r. 
We see that w; has internal layers. In Figure we plot the graph x — ► t> (r, x)/r for r = 60. We 
first see that the function takes all its values in a small interval and has very rapid variations 
with respect to x (is nearly discontinuous). This does not contradict the theory, because there 
are no uniform estimates on the modulus of continuity of v(t, -)/t. 
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Figure 3. First case, p = 1.3: w(t, 0,0) - (w(t)) (left) and v(t, 0) - {v(r)) 
(right) function of r 



4.1.2. Second case. We consider a two dimensional problem, where the Hamiltonian is 
H(x,u,p) = cos(27rxi)+cos(27r2;2)-|-cos(27r(xi— X2))+sin(27rti)+ 



COs(27TXl) sin(27TX2) 

2 4 



\p\- 



For this case, only the Imbert-Monneau cell problems have been approximated on uniform grids 
with step 1/200. The time step is 0.005. In Figured! we plot the contours and the graph of 
the effective Hamiltonian computed with the high order method. We can see that the effective 
Hamiltonian is symmetric with respect to p = (0,0), constant for small vectors p. In Figure [7J 
we plot as a function of r. We see that this function converges when r — > oo. In Figure 

El we plot the contours of v(t, -)/t for r = 59.935 and p = (1, 1). We see that for large values of 
r, v is close to discontinuous. 



Appendix A. 



Proof of Lemma 12.51 To show that the sequence is convergent it suffices to show that for any 
s E R ipe' S (s) is a Cauchy sequence. Fix s el and let iq € Z be the closest integer to s, i.e., 
s = ioe + 7e, with 7 6 (— |, ^] . Let k > m > |io|, then, by assumptions (|2.8p we have 



-m— 1 



E 

i=—k 
—m—1 

E 

i=—k 



s — ei 



s — ei 



+ E 



i=m+l 



s — ei 



e(k — m) 



+ E 



-m— 1 



i=m+l 
fc 



s — ez 



£^ 2 E^ + ^ 2 E 77^ 



i=—k 
sr2 



Similarly, it can be showed that 



i=m+l 
s;2 k 



S — ei) 



c (i - » + 7) 2 



+ *T E 



m,<5 



» > E 



e . ' (i - « + 7) 5 

i=m+l 



x2 1 
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"w(60, . , . ) /60" 




"w(60,0.13, .)/60" 




drri Brr. 



-e-. 25 



Figure 4. First case, Barles cell problem, p = 1.3. Top: contour lines of 
w(t, -)/t on a period as a function of (x,y). Bottom: the cross-section x = 0.13. 



Hence \ife' S (s) — ipT'^is)] — > as m,i — > +oo. Similar arguments show that the sequence 
(<pi' n )' converge uniformly on compact sets of R. This implies that <pf is of class C 1 with 
(^)'( S )=lim n _ >+00 (^)'( S ). 

Now, let us show (|2.9jl . Let s = io^ + 7^ for some G 2 and 7 G [0, 1). Then 
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Figure 5. First case, Imbert-Monneau cell problem, p = 1.3: Top: Third order 
Runge Kutta/WENO scheme: v(t,x)/t as a function of x for r = 60; the right 
part is a zoom. Bottom: same computation with Euler/Godunov scheme with 
the same grid parameters: some oscillations are smeared out, but the average 
value of the solution is well computed. 
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Figure 6. Second case, the effective Hamiltonian computed by solving Imbert- 
Monneau cell problems. 
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Figure 7. Second case, p = (1, 1). The median value of v (r, -)/r on a period as 
a function of r. 



Similarly 



(5 2 n+i ° I 



-*oe>e[<*(-^)-l] - 7 ^E 



j2 , 

T 2 2^ 



- (i + 7 )2 



- (i - 7 )2 



Letting n — > +oo, we get 

tp 5 e(s) -ioe-e 



2 +°° 



- (i - 7)2 
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Figure 8. Second case, the contours of the solution of Imbert-Monneau cell 
problem for p = (1, 1) at time r = 59.935. 

If 7 > then <p (^) - 1 -> as 5 -» 0+ and p£(«) -> i e if 6 -> 0+. If 7 = 0, then 
~» (*o - l)e + </»(0)e if 5 -» 0+ and ([IS]) is proved. □ 
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